Gut microbiome reflect adaptation of earthworms to cave and surface environments

Background Caves are special natural laboratories for most biota and the cave communities are unique. Establishing population in cave is accompanied with modifications in adaptability for most animals. To date, little is known about the survival mechanisms of soil animals in cave environments, albeit they play vital roles in most terrestrial ecosystems. Here, we investigated whether and how gut microbes would contribute to the adaptation of earthworms by comparing the gut microbiome of two earthworm species from the surface and caves. Results Two dominant earthworm species inhabited caves, i.e., Allolobophora chlorotica and Aporrectodea rosea. Compared with the counterparts on the surface, A. rosea significantly decreased population in the cave, while A. chlorotica didn’t change. Microbial taxonomic and phylogenetic diversities between the earthworm gut and soil environment were asynchronic with functional diversity, with functional gene diversity been always higher in earthworm gut than in soil, but species richness and phylogenetic diversity lower. In addition, earthworm gut microbiome were characterized by higher rrn operon numbers and lower network complexity than soil microbiota. Conclusions Different fitness of the two earthworm species in cave is likely to coincide with gut microbiota, suggesting interactions between host and gut microbiome are essential for soil animals in adapting to new environments. The functional gene diversity provided by gut microbiome is more important than taxonomic or phylogenetic diversity in regulating host adaptability. A stable and high-efficient gut microbiome, including microbiota and metabolism genes, encoded potential functions required by the animal hosts during the processes of adapting to and establishing in the cave environments. Our study also demonstrates how the applications of microbial functional traits analysis may advance our understanding of animal-microbe interactions that may aid animals to survive in extreme ecosystems. Supplementary Information The online version contains supplementary material available at 10.1186/s42523-022-00200-0.


Background
Terrestrial caves differ from surface habitats and are regarded as "natural laboratories" [1]. Organisms in caves are subjected to strong selective pressures that are rather different from surface, such as constant dark and static climatic conditions [2]. Further, food resources for animals are depleted in cave as compared to the surface, which rely on photosynthesis for primary production [3]. Instead, food webs in caves are based on microbes, making the interactions in soil food webs in cave hitherto unknown [4]. Cave-dwelling animals therefore have to modify their feeding strategies in ways different from surface animals to cope with food deficiency [5].
On the surface, earthworms are common ecosystem engineers which provide a variety of ecosystem functions and services [6]. Some of them are able to establish stable populations in subterrain caves [7,8]. Species that survive in caves need to adjust to the rather different environment compared with those living in other ecosystems. Compared to cave arthropods that are usually pale and blind [9], cave earthworms appear to be similar in external morphological characters to their relatives inhabiting surface layers. However, adaptive strategies of earthworms to cave environments and the manner how different earthworm species sustain populations in the resource-limited environments remain obscure.
Gut-associated microbes are considered as "nutrient factories" that increase host fitness in many ways [10][11][12]. With the help of diverse functional genes encoded by gut microbes, the hosts may digest a wide range of compounds, thereby surviving in unfavorable environments [13]. Reciprocally, gut microbiome are regulated by the physiological conditions and feeding diets of the hosts [14,15]. For example, the community structure and functional potential of the gut microbiome differ in gut compartments and are associated with the feeding strategies of the hosts [16]. Besides, species-specific effects from the hosts on gut microbiome have also been discovered for soil animals [17]. Therefore, cave earthworms are presumably to harbor specific gut microbiome which may help them to adapt to the cave environment, and the same earthworm species may reshape their gut microbiome when living in a different environment.
Here, we asked whether gut microbiome of cave earthworms differ from the respective earthworm species inhabiting surface and if the gut microbiome may provide functions for earthworm adaptation to the optimal environments (i.e., surface or cave). We explored patterns of gut bacterial communities and predicted functional genes encoded by gut bacteria as well as the copy numbers of the 16S rRNA gene to reflect the ecological strategies of bacterial communities in nutrient exploitation [18]. We also examined networks of gut bacterial communities which suggest interactions and stability of the gut microbiome [19]. We hypothesized that 1) gut bacterial communities differ between cave and surface earthworm populations in taxonomic, functional and phylogenetic diversity and the difference depends on the species adaptability to environments; 2) The assembly processes of gut microbiome are more deterministic than that of soil microbiota since cave earthworms are likely to select bacteria of certain functions for the adaptation in the specific environment; 3) Gut microbiome are characterized by less fast growing species and more stable networks when the earthworm hosts inhabit favorable environments.

Study sites and sampling
The study sites are located in two interconnected cave systems, the Amatérská Cave and the Sloupsko-Šošůvské Caves, in the Moravian Karst Protected Landscape Area in the south-east part of the Czech Republic ( Fig. 1 and Table 1). The cave system is associated with streams, Sloupský potok and Bílá voda, which later merge to the Punkva River. The water mainly flows underground, dropping to the lime bedrock massif. The gallery-like caves were originally formed by streams and have been connected with the surface layers via water flow. Nowadays they are situated in the vadose zone far from the streams, with the only water supply being infiltration of rainwater through soil and bedrock rifts, except for extreme flooding events. The two sampled cave systems are separated by a series of water siphons and are accessible only through an artificial corridor. No bats are present here and support of organic matter is possible via floods only. Just the corridors of the Sloupsko-Šošůvské Caves (Site 1; Fig. 1 and Table 1) are connected directly with the surface. The sampling took place in the center of the caves, which were far away from the water siphons. Soils in all sites are of grey rendzina type, affected by close permanent or semi-permanent water stream. Vegetation of surface sites is very similar, represented by Stellario nemorum-Alnetum glutinosae. Soil pH values of the cave and soil substrates ranged between 8.5 and 7.0, soil organic carbon (SOC) 7-17.1 g/kg. [20].
Earthworms and soils from each site were sampled during the data 3-10 May 2016, which was spring for the sampling sites and corresponding to the growing season for the earthworms, evidenced by the casts they produced in the field. The sampling sites were surveyed for common soil animals, e.g., mesofauna springtails and macrofauna earthworms, and the present study was focused on the earworms. In the cave, the soil was mixed with earthworm casts (Fig. 1A), thus the sampling of cave soil was a mixture of soil and earthworm casts. On the surface, the soil sampling was conducted in an area dominated by a specific earthworm species. When sampling, five 5 × 5 m quadrats were randomly set up at each site with a distance of at least 10 m (Fig. 1B, C), and eight soil cores were collected and mixed for one composite sample for each quadrat. The soil cores were taken with a 2.5 cm diameter cylinder, and to a depth of 10 cm, or until reaching the rock. This method enabled at least 100 g of soil for each core. After soils were collected, the earthworms were dug out from each quadrat. In each quadrat, earthworms were hand-sorted and preserved in their living soils and were transported with ice to the lab. In the lab, the earthworms were identified to species level and stored in absolute ethanol prior to molecular gut content analysis. The abundance of the earthworms per site was calculated as the mean of earthworm density in the eight quadrats. After sampling, soils with earthworms submerged in their living soils were transported with ice to the lab. In the lab, the earthworms were fixed and stored in absolute ethanol and then identified to species level prior to molecular gut content analysis. Only two earthworm species were found in the cave systems, i.e., Allolobophora chlorotica and Aporrectodea rosea in Sloupsko-Šošůvské Cave (Fig. 1, site1) and Amatérská Cave ( Fig. 1, site 2), respectively. A. chlorotica and A. rosea belong to the same family, i.e., Lumbricidae, and are both widespread in Europe [21].

Molecular gut content analysis
Five individuals of earthworms from each site were dissected and separately used for DNA extraction during the molecular analysis. The earthworms were dissected aseptically under a stereomicroscope. An incision was made longitudinally along the body wall and the whole gut, from the clitellum to anus, was removed and placed in a 1.5 mL Eppendorf tube. Thereafter, total DNA of the gut content as well as soils were extracted using the FastDNA Spin Kit for Soil and the FastPrep Instrument (MP Biomedicals, Santa Ana, CA, USA). All steps were carried out following the manufacturer's instructions. The quality and quantity of the extracted DNA were certified with 1% agarose gel electrophoresis and a Nanodrop-2000 spectrophotometer (NanoDrop Technologies Inc. Wilmington, DE, USA), respectively. The V4 hypervariable region of the bacterial 16S rRNA gene was amplified and sequenced with a Miseq sequencer at the University of Illinois-Chicago with primers of 16S rRNA gene V4 region (FWD: 5'-GTG YCA GCMGCC GCG GTAA-3'; REV: 5'-GGA CTA CNVGGG TWT CTAAT-3') were used following the EMP protocol (https:// earth micro biome. org/ proto cols-and-stand ards/ 16s/). Negative controls that replace DNA templates with sterilized water were included in the amplification period. The raw sequences were deposited in NCBI Sequence Read Archive under the accession number PRJNA400302.

Sequence data processing
Paired-end sequence data were joined, demultiplexed and analyzed using the QIIME 1.9.1 pipeline [22]. Briefly, sequence lengths < 200 bp, of average quality score < 20 or with ambiguous characters were discarded. After chimeras and singletons were removed, closed reference operational taxonomic units (OTUs) were clustered on the basis of 97% similarity. Taxonomy of bacterial OTUs was assigned using Greengenes v13_8. A phylogenetic tree was generated using "make_phylogeny.py" by the default setting of the "FastTree" method. The resulting Table 1 Abundance of earthworms (Allolobophora chlorotica and Apporrectodea rosea) and the soil properties in caves and surfaces. Values are presented as mean ± SD ┼ Locations of sampling sites as shown in Fig. 1

Habitat
Cave Surface OTU table was then rarefied to 9800 sequences per sample before further analysis.

Statistical analysis
The 16S rRNA gene copy numbers as well as functional gene abundance were calibrated and predicted by PIC-RUSt [23]. The abundance of functional genes was predicted using a script predict_metagenomes.py implemented in the PICRUSt according to the recommended protocol. The predicted genes were then grouped at the first KEGG level using the script categorize_by_function. py implemented in the PICRUSt. Other statistical analyses were performed in R 4.0.0 [24]. Student's t-test was used for comparing the mean abundance of earthworms between cave and surface. The standardized effect size of abundance weighted mean phylogenetic distance of the bacterial community was quantified using the function ses.mpd implemented in the R package "picante" [25]. OTU numbers, diversities and phylogenetic relatedness of microbial communities were compared between treatments using ANOVA followed by Tukey's HSD test. Community weighted means (CWM) of the 16S rRNA gene copy numbers were calculated with the equation: CWMoperon = n 1 Pi * mi , where Pi and mi is the proportion and operon numbers of each bacterial OTU. Pairwise correlations of the bacterial OTUs within treatments were calculated using the command sparcc with 1000 bootstraps in the program mother v.1.35.0 [26]. Significant correlations were set by R 2 > 0.7 and P < 0.01. Topological properties of the bacterial community network of each treatment included (I) numbers of nodes and edges, (II) average degree, which measures network complexity, and (III) average path length (i.e., distance between any two nodes). Network properties were calculated using the "igraph" package in R [27].

Microbial diversity
Functional diversity inferred by the predicted functional gene richness was greater in the gut of earthworms than those in the soils (Fig. 2A). Gut bacteria of both A. chlorotica and A. rosea were functionally more diverse in the caves than in the surface ( Fig. 2A). Regarding the eight categories of predicted functions, both A. chlorotica and A. rosea exhibited greater genes related to metabolisms (Additional file 1: Figure S1). Gut microbiome holds more abundance of functional genes than that of soils, but the difference of genes between cave and surface was not significant (Fig. 2A). Taxonomic and phylogenetic diversities of bacterial communities, however, were lower in the gut of earthworms than those in the soils (Fig. 2B and C). Regardless of gut or soil microbiota, the taxonomic and phylogenetic diversity was lower in the caves than in the surface for both A. chlorotica and A. rosea ( Fig. 2B and C), with the greater reductions for A. rosea (F = 4.93 and 14.77 for A. chlorotica and A. rosea, respectively; P < 0.05).

Phylogenetic relatedness
Except for the bacterial communities in the gut of A. chlorotica from the surface, which exhibited a random pattern of phylogenetic relatedness, the bacterial communities of all the other treatments showed phylogenetic clustering (Fig. 3). The standardized effect size of mean phylogenetic distance of soil bacteria was significantly lower than in the gut of A. chlorotica irrespective of the habitats (i.e., surface or cave; P < 0.05). However, in the gut of A. rosea bacterial communities were more phylogenetic clustered than the soil bacterial communities if they were collected from the surface but not the cave. For A. chlorotica, the mean phylogenetic distance of the gut bacterial communities was greater in the surface than cave, while for A. rosea, the reverse was true (Fig. 3B).

Community weighted mean operon numbers
For both earthworm species in general, the communityweighted mean of the 16S rRNA gene (rrn) copy numbers of bacteria were significantly higher in the gut of earthworms than that in the soils (mean value 3.8 and 2.5, respectively; Fig. 4). The CWM operon numbers were not different between the surface and cave soils (P > 0.05). However, in the gut of A. chlorotica, the value was significantly greater in the surface than in the cave, while an opposite pattern was found in the gut of A. rosea (P < 0.05).

Co-occurrence networks of microbial communities
Networks mainly consisted of the most abundant phyla and comprised highly connected OTUs structuring densely connected groups of nodes ( Fig. 5 and Table 2). The networks of bacterial communities in the gut of both earthworm species exhibited fewer degrees. The degree of the network in the gut of A. rosea was reduced ~ 50%, while A. chlorotica increased 18% when living in the cave as compared to the surface layers. Gong et al. Animal Microbiome (2022) 4:47 Discussion A. chlorotica and A. rosea exhibited different variations of abundance in caves. Presumably, they might consume more forms of food with the aid of their gut microbes. Adaptation of animals to different environments usually requires physiological adjustments, including changes in biochemical activities. However, physiological or geneticbased adaptation of the earthworms usually takes generations, longer than the lifespan of individual earthworms. The gut microbiome, as part of the holobionts, may facilitate the animal hosts to adapt to unfavorable environments by a diverse encoding of genes [28,29]. The gut microbiome of earthworms, therefore, could instantly be recruited by the host and help them to utilize more diverse food compounds, thereby increasing their fitness in different environments [30].

Functional diversity of microbiota in the earthworm gut
We found that cave earthworms harbor more diverse functional genes in their gut, despite lower taxonomic and phylogenetic diversities, supporting our first hypothesis. Both earthworms enriched the genes related to the function of metabolism in their gut microbiomes, which is likely to provide essential metabolites to them and increase their survival rates [31]. Higher proportions of metabolism-related genes were found in the gut of soil animal microbiota compared to other functions.
It is evidenced soil animals are highly dependent on the metabolism genes from their gut microbiome for nutrient needs [32,33]. As the cave environment is deficient in food sources, when the earthworms are living in the caves, their gut might be stimulated to serve as a more efficient "nutrient factory" [12]. Notably, a large proportion of the genes, especially from the gut microbiome of earthworms in the caves, was unclassified. This demonstrates that the gut microbiome of cave earthworms is more functionally diverse than has been seen. More studies integrating new technologies should be conducted to uncover their roles.

Deterministic community assembly in the earthworm gut
Our results show that the assembly processes of earthworm gut microbiome were deterministic, supporting our second hypothesis. The oxygen and water contents differed along the digestive tract of earthworms, thus the profile of microbiota is shaped by the digestive environment [34][35][36]. In addition, studies have revealed that the food source is another deterministic factor shaping the gut microbiome [37][38][39][40], both phylogeny and food preference of the host may deterministically shape the gut microbiome of soil animals [41]. For surviving in caves, earthworms therefore might select microbiota coding more metabolic activities.

Community features of gut microbiome
The fact that gut microbiome holds greater numbers of CWM copy numbers of the 16S rRNA gene than soil microbiota suggesting a more nutrient demand in the gut microbiome. The multiplicity of rRNA genes is an indicator of the ecological strategy of bacteria for nutrient exploitation [18]. For example, during the exponential growth phase, the number of rRNA operons of Escherichia coli may increase from 7 to 36 [42,43]. Therefore, communities dominated by bacteria with fewer rrn copy numbers usually have a higher nutrient use efficiency than high-rrn-dominated communities [44]. In caves, where soil organic matter is different compared with surface systems, earthworms might benefit from the help of a highly efficient gut microbiome supporting higher metabolism [45]. That was the reason A. chlorotica established a more stable population, compared with A. rosea.

Conclusions
The present study demonstrates the functional roles of gut microbiome in contributing to host adaptations of two earthworms in surface and cave environments. Our results reflect the tight interactions between host earthworms and their gut microbiome. The gut microbiome exhibits a more functional diversity in caves, which can be interpreted as evidence of stronger food limitation. A more stable and highly efficient microbiota providing metabolites is needed for the earthworms to survive in the resource-limited cave habitat. Together, the gut microbiome-host crosstalk is of pivotal importance in facilitating the animal hosts in their physiological adaptation and even the population expansion. Copy numbers were estimated using PICRUSt and weighted values were obtained by multiplying copy numbers by the relative abundance for each operational taxonomic unit and taking the sum of these values for each community